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I INTRODUCTION 
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The  accuracy  and  reliability  of  data  from  microwave  radars 
can  be  improved  if  amplitude  and  angle  corrections  are  applied 
to  the  received  radar  signals.  Such  corrections  can  be  computed 
if  the  radio  refractive  index  profile  along  the  various  propa- 
gation paths  is  available. 

Naval  studies  have  shown  that  the  most  difficult  problem  in 
reliable  mapping  of  the  atmospheric  radio  refractive  index  is  the 
forecasting  of  the  water-vapor  distribution  along  the  line  of 
sight  between  the  target  and  the  transmitter.  Accurate  forecas- 
ting of  the  water-vapor  concentration  is  difficult  because  of 
inadequate  knowledge  of  the  water-vapor's  distribution  at  any 
given  time  and  because  of  the  difficulty  in  predicting  vertical 
atmospheric  motions  which,  to  a large  extent,  control  the  changes 
in  the  water-vapor  distribution.  The  verical  motions  that  cause 
major  changes  in  the  water-vapor  distribution  are  usually  small, 
cannot  generally  be  measured  and  must  be  inferred  from  large  scale 
horizontal  flow  fields. 

The  simultaneous  measurement  of  temperature  has  been  identi- 
fied as  a significant  secondary  parameter  of  interest  to  the  Navy. 
Present  techniques  for  the  measurement  of  water-vapor  and  temper- 
ature in  the  atmosphere  utilize  radiosondes  (ballon  launched  and 
air  dropped) . In  addition  to  requiring  the  use  of  expendables, 
radiosonde  techniques  are  deficient  because  they  are  volume  and 
direction  limited.  The  volume  of  atmospheric  coverage  of  radio- 
sondes is  not  sufficient  to  perform  adequate  mapping  of  spatial 
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and  temporal  features  important  to  radio  meteorology.  The  radio- 
sonde vertical  profiles  are  too  few  and  far  between  and  do  not 
sample  in  real  time  the  volume  of  the  atmosphere  for  which  re- 
fractive index  corrections  are  required. 

Although  currently,  the  Navy  is  investigating  the  use  of  Lidar 
as  a viable  means  of  determining  the  vertical  water-vapor  dis- 
tribution in  the  atmosphere,  this  study  is  addressed  to  the  ab- 
ility of  passive  infrared  and  microwave  radiometry  to  extract 
the  same  information. 

The  relative  merits  of  active  and  passive  methods  are  not  the 
sub'  of  this  study;  however,  it  is  believed  that  ultimately 
t ' ve  system  will  be  found  to  be  less  expensive  and  bulky 

he  Lidar  system  while  providing  adequate  mapping  of  the 

I 

water-vapor  distribution  necessary  for  radar  propagation  correc- 
I tions. 

In  order  to  assess  the  adequacy  of  the  passive  method,  section 
i * 

j II  analyzes  the  degree  of  spatial  resolution  required  to  adequate- 

ly describe  radio  wave  propagation. 

1 
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II  SPATIAL  RESOLUTION  REQUIREMENTS  FOR  DESCRIBING  RADIO  WAVE 


■ PROPAGATION 

This  section  establishes  the  guidelines  for  determining  the 
spatial  resolution  requirements  necessary  to  adequately  describe 
radio  wave  propagation  through  the  atmosphere  whose  water-vapor 
and  temperature  distributions  are  changing  in  time  and  space. 
These  guidelines  are  necessary  to  properly  assess  whether  a pas- 
sive system  can  adequately  recover  the  water-vapor  and  tempera- 
ture profiles.  The  degree  of  recovery  required  will  be  dictated 
by  the  degree  of  recovery  necessary  to  adequately  describe  the 
radio  wave  trajectory. 

Due  to  the  variety  and  magnitude  of  meteorological  variables 
involved,  we  have  made  the  following  assumptions  which  are  con- 
sidered realistic; 

(a)  All  refractive  index  gradients  are  normal  to  the  earth's 
surface  and  are  a function  of  the  radial  distance  only 

(b)  The  verical  gradients  do  not  change  with  horizontal 
extent  (for  distances  to  10  Km;  this  is  considered 
reasonable  for  a marine  environment) 

(c)  Radio  wavelengths  of  6mm  or  greater  are  being  utilized 

(d)  The  radio  wave  refractive  index  is  affected  by  tempera- 
ture, water-vapor  and  pressure  only  (i.e.  cloudless) 

(e)  The  radio  wave  propagation  can  be  described  accurately 
by  geometric  optics 
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with  these  assumptions  in  mind,  we  proceed  to  derive  the  radius  ’ 

V ' 

of  curvature  of  a radio  wave  ray.  ; 

A.  Derivation  Of  The  Radius  Of  Curvature 

The  vector  differential  form  of  the  path  of  rays  in  an  inhomo- 

(1) 

geneous  medium  is  given  by, 

4^)  - Vn  (1) 

<JS  ' 

where  s=  distance  along  ray 

n=  refractive  index,  a function  of  r 
r=  radius  vector 
Vh=  gradient  of  n 
Expanding  we  have  that, 

. ilT  n ^ 

ds  ds  ds^ 

(2) 

/\ 

We  identify  dr/ds=  t(unit  vector  tangent  to  the  s direction, 
see  Figure  1) . 

Consequently, 

^ 4%-  ^ ^ 

For  a spherically  symmetrical  medium  in  which  V'n  = r'  dn/dr 
(r  the  unit  vector  in  the  r direction) , assumption  (a)  above, 
we  have  that, 

t In.  4-  ^ (4) 

d s d s*"  cl  r 

4 


I 


We  identify  the  term  with  iT,  the  curvature  vector,  whose 

magnitude  is  1/j^  , Ji  the  radius  of  curvature  and  whose  direction 
is  perpendicular  to  't  . 

After  rearranging,  equation  (4)  becomes, 

- dn  — d*^ 

15) 

Multiplying  through  by  iT  we  have, 

(6) 

As  is  perpendicular  to this  becomes, 

(7) 


n _L.  =.  dn L Cos ^ 

6t-  p 


6t-  jo 

with = the  angle  between  'r  and  K.  Since  n is  close  to  one. 


we  have  that. 


djl_  coS  A 


From  Figure  1,  oi  / where  ^ is  the  angle  the  ray  makes 
with  the  local  horizon.  Finally  we  have  that, 

^ cos  0<.  ,1, 

For  a spherically  symmetric  distribution,  the  ray  is  a part 
of  a circle  whose  radius  is and  is  confined  to  a plane  formed 
by  the  initial  ray  direction.  The  derivation  of  the  coordinates 
of  the  center  of  the -circle  is  found  to  be, 

0-  — S»r»  0(.  ^2] 

b - x-- yOCoSoC 
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Consequently,  given  0<.  and  dn/dr  (constant  gradient)  with  the 
earth's  center  at  the  center  of  the  coordinate  system,  the  ray 
trajectory  follows  the -circle  whose  center  is  at  +a  and  -b 
for  a negative  gradient. 

B.  Derivation  of  the  Angular  Variation  of  the  Ray 

f^ct  that  the  ray  is  confined  to  move  along  the  arc  of 
a great  circle  allows  the  derivation  of  the  relationship  between 
T /the  angle  the  ray  makes  with  the  local  horizon,  and  the  known 
variables  using  only  geometric  considerations.  Referring  to  Figure 
1,  using  the  law  of  sines  we  have  that. 


X = “K  6 _ 

Sin  S 

where  6=  180- €.  with  €=  cX  + 0 . Consequently, 

X ■=  "R 

Sin  € 

Similarly, 


(13) 


(14) 


Sin  r = sin  e 

Substituting  from  equation  (14)  we  have  that. 


(15) 


Sin  r = ^ 


Sin  C«  + 6)  — T^SinG 


(16: 


Expanding  and  substituting  the  expression  for  derived  in 
section  lA,  we  obtain. 


S(nr  = sin0cos<xCl  s«no(  cos  G 

an  (17 ) 


As  we  are  concerned  with  horizontal  distances  only  up  to  lOKM, 
then  Sin  G=  © and  CoS  0 1 resulting  in. 
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Yc  0 cos  o<  ( ' ^ ^ 


(18) 


Sin 


Equation  (18)  implies  a negative  gradient  and  will  be  used  below 
to  calculate  the  spatial  resolution  requirements. 

C.  Derivation  of  Spatial  Resolution  Equations 
Before  equation  (18)  can  be  used  to  determine  the  spatial 
resolution  requirements,  expressions  for  dn/dr  must  be  derived 
in  terms  of  variations  in  pressure,  temperature  and  water  vapor 
pressure . 

The  instantaneous  refractive  index  can  be  written  as, 

n = Ho  -b  S (19) 

where  n^  = the  refractive  index  of  some  reference  set  of  con- 
ditions (e.g.  no  atmosphere  or  a standard  atmosphere  with  60% 
relative  humidity)  and  the  variation  from  the  reference. 

Letting  dn  “^An  we  have  that, 

+ <S  (n^.o)  ~ n,,o  - S(ni,o)  (20 ) 

or, 

AnsAHo-*-  (21) 


However, An^  can  be  expressed  in  terms  of  the  refractive  index 
modulus  N as  used  by  meteorologists  where, 

An^  =■  ANo  I 0“^  ( 


Consequently , 

6(Ano)=  ^ CaNo)  k 

Substituting  in  equation  (21)  results  in. 
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n={ 


S (^No)| 


X lO 


-fe 


(24) 


Additionally,  as  dr'^^H,  the  change  in  height,  we  have  that, 

^ An  - •*“  ^ lo"^ 

"dT  A.H  <25 

We  call  ANq/  A HxlO  K , as  it  represents  a constant  gradient 
or  no  gradient.  Consequently,  R dn/dr  becomes, 

R-dO- RK  + 4.3T1 

d>-  (26 

( 2 ) 

where  we  use  the  average  earth  radius  as  6.371x10^  meters 
Substituting  into  equation  (18)  and  expanding  we  obtain. 


SinTr  ecoSo(  -eRKCoS«-b^52l.ecoSO(£C^^^i^^  ^27) 

A H 

As  we  are  interested  in  the  deviation  of  T f rom  a reference  grad- 
ient we  obtain  that, 

A (smr)  s«nX  -sinr^=  6:22i©i2S2Li^\28) 

A H 


where , 

"Yo  = ray  angle  with  reference  fixed  gradient 

angle  with  fixed  and  variable  gradient 
Equation  (28)  relates  the  deviation  in  Y*  due  to  a deviation  in 
N from  the  fixed  reference.  In  order  to  put  equation  (28)  in  a 
more  useful  form,  we  note  that  a cosO(  = H cosO</R  tano^  of  the 
disturbance  of  extent  H.  Substituting  we  have  that. 


A (smr-)  = H &CAN.)co.5^  ^ , Q-fc 

^ w +an  oc 


(29) 


We  will  call  the  expression  H <S(ANq)/aH,  t*  , a measure  of 
the  intensity  of  the  disturbance  causing  the  deviation  from  the 
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reference.  As  we  are  dealing  with  small  changes  in  angle,  we  equate 
^ (sinT)  to  d(sinY')  which  equals  cosVdT"  or  cos  VAV*. 

Substi- 
tuting in  equation  (29)  and  noting  that  cos  0^  ~ cost*  we  finally 
have  that. 


This  relates  the  strength  and  extent  of  the  disturbance  to  the 
change  in  that  can  be  tolerated  and  the  angle  of  elevation  of  the 
ray.  Before  equation  (30)  can  be  used  to  determine  the  spatial  res- 
olution requirements,  an  estimate  must  be  made  of  . 

D.  Effect  of  Meteorological  Variables  on  N 

Under  assumption  (c)  of  part  A,  the  refractive  index  modulus 
(3) 


can  be  written  as 


- 77.  6 P 573. 25-6  "Pw 

-r 


where , 

P = atmospheric  pressure  in  millibars 
T = temperature  in  °K 
P^=  water  vapor  pressure  in  millibars 
As  we  are  interested  in  the  AN  due  to  these  variables  then. 


AP  +-  APw  (32) 

dP  37^  O vw 

Carrying  out  the  differentiation,  we  obtain, 

/V  Kt  = 77-6  AP  _ jnTfeP  ^74fe.5ia.PwlAT  4-  (33) 

I J T' 

For  heights  up  to  2 KM,  the  variation  in  P,T,  and  P are  app- 

w 

roximately , 


9 


1020<  P>  750 


315<  T > 270  (34) 

30<P>  3 

w 

For  hot  humid  air, 

^Mr.25A'P  4“  3.76  ATJ,/  — l-SIAT  (35) 

and  for  cool  dry  air, 

-hS.TI  AT^  —.91  AT  (36) 

As  we  are  interested  in  estimates  of  the  spatial  resolution 
requirements,  we  average  the  two  extremes  and  let  the  devia- 
tion be  represented  by, 

<A^i>  =--27  AT  ■m.44APU/  ^37) 

This  infers  that  the  AN  change  with  P,T  and  P^  is  constant  up 
to  2 KM,  As  the  effect  of  pressure  relative  to  T and  P^  is  small, 
for  the  estimation  of  the  spatial  resolution  we  will  use  the  ex- 
pression, 

— 4-^4APw  ■“  AT 

Dividing  this  by  the  A H in  which  this  change  occurs  and  mul- 
tiplying by  H the  physical  extent  of  the  disturbance  provides  the 
^ to  be  used  in  equation  (30). 

As  an  example  of  the  calculation  of  the  magnitude  of  we 

consider  an  inversion  layer.  An  inversion  layer  usually  indicates 
the  presence  of  a humid  mass  under  a dry  one.  The  transition  from 
humid  air  to  dry  air  air  causes  the  largest  change  in  N.  For  ex- 
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ample,  a temperature  inversion  of  2°/100  meters  located  near  the 
ground  would  produce  a change  of  = -.3  and  a T = .02. 

From  equation  (38),  = -1.  33  - .02  = -1.35/meter.  If  the  phy- 

sical extent  of  the  disturbance  is  200  meters  (typical  of  the  Med- 
iterranean area)  then  ^ = -270.  The  minus  sign  indicating  the 
direction  of  change  from  the  reference.  If  more  than  one  distur- 
bance is  present,  the  the  effect  is  additive  as  the  displacement 
effect  is  negligible.  Consequently,  we  consider  a range  of  = 1 to 
10^.  = 10^  representing  very  unusuall  conditions  and  must 

be  considered  atypical.  If  an  acceptable  A'Y'is  known,  the  spatial 
resolution  can  be  determined  from  equation  (30). 

E.  Calculation  of  the  Spatial  Resolution  Requirements 

In  Figure  2,  we  have  plotted  log  t against  0(  for  three  angles 
of  Air  considered  appropriate  for  radio  wave  propagation.  These 
curves  were  calculated  from  equation  (30). 

We  consider  first  what  spatial  resolution  is  required  to  adequat- 
ely determine  a standard  atmosphere  with  a 60%  relative  humidity. 
Calculation  indicates  that  the  A n/  A H = -.04/meter.  If  we 
use  as  a reference  no  gradient  (straight  line  trajectories) , we 
can  treat  the  standard  atmosphere  gradient  as  the  disturbance  with 
a physical  extent  of  2 KM.  Consequently,  ^ = 80.  Referring  to  Fig- 
ure 2,  if  the  tracking  error  is  to  be  .01°  or  less,  we  find  that 
for  an  elevation  angle  of  25°,  the  presence  of  the  atmosphere  does 
not  effect  the  straight  line  trajectory  to  within  .01°.  For  eleva- 
tion angles  less  than  25°,  the  spatial  resolution  required  will  be 
determined  by  . For  example,  at  0^  = 1°,  ?“  = 3.  Hence  in  order 

to  predict  the  trajectory  to  within  .01°,  spatial  measurements  of 
about  2000/80/3  or  about  every  74  meters  is  necessary. 
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For  inversion  layers  and  humid  air  masses  where could  app- 
roach 1000  or  more.  Figure  2 indicates  that  for  a .01°error  and 
a 45°  elevation  angle,  the  spatial  resolution  required  is  1000/ 
1000/175=  175  mete..  if  the  disturbance  has  a physical  extent  of 
1000  meters.  Other  vation  angles  require  different  spatial  re- 
solution measurementij  to  define  the  gradient. 

F.  Analysis 

The  examples  cited  above  indicate  that  a single  spatial 
resolution  does  not  suffice  but  depends  on  the  angle  of  elevation, 
the  gradient  caused  by  the  disturbance  and  the  error  in  refracted 
angle  that  can  be  tolerated. 

Additionally,  for  elevation  angles  30°  or  greater  with  a 
tolerable  angle  of  .1°  or  greater;  no  set  of  inversion  requires 
any  better  spatial  resolution  measurements  than  are  necessary  for 
the  standard  atmospheric  gradient. 

Unless  trajectories  close  to  the  horizon  are  required,  we 
conclude  that  coarse  spatial  resolution  measurements  are  all  that 
are  necessary  to  adequately  describe  the  radio  wave  propagation 
for  most  applications. 
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RAY  TRAJECTORY 


FIGURE  1.  GEOMETRY  USED  TO  DERIVE  RAY  EQUATIONS.  (SEE  TEXT) 


Ill  MICROWAVE  AND  INFRARED  TRANSFER  EQUATIONS 


In  order  to  arrive  at  an  adequate  method  for  recovering  both 
the  temperature  and  water-vapor  profiles,  we  examine  the  atmos- 
pheric radiation  due  to  water-vapor  in  the  infrared  and  micro- 
wave  regions  and  due  to  atmospheric  oxygen  in  the  microwave  re- 
gion and  carbon  dioxide  in  the  infrared  r;igion.  Although  the  car- 
bon dioxide  radiation  could  be  used  to  determine  the  temperature 
profile,  it  is  not  considered  in  this  study  as  it  is  believed  that 
the  infrared  sensors  at  the  carbon  dioxide  radi,ation  wavelength 
(above  15  microns)  are  not  suitable  as  those  at  the  water-vapor 
radiation  wavelength  (6  to  7 microns) . In  addition,  as  the  micro- 
wave  technology  is  so  highly  developed,  it  would  appear  that  it 
is  a more  suitable  region  of  the  spectrum  to  use  in  developing  a 
sensor  for  temperature  and  water-vapor.  As  both  carbon  dioxide 
and  oxygen  can  be  considered  to  have  a constant  mixing  ratio  up 
to  2 Km.,  either  could  be  used  to  infer  the  temperature.  On  bal- 
ance, for  this  study,  we  have  chosen  to  utilize  the  oxygen  rad- 
iation for  the  temperature  profile  determination. 

A.  The  Radiative  Transfer  Equation 

The  amount  of  radiation  emitted  by  the  atmosphere  and  received 
by  a ground-based  radiometer  can  be  determined  by  the  intergration 
of  the  transfer  equation  along  the  ray  path.  If  we  consider  the 
atmosphere  to  be  in  local  thermodynamic  equlibrium  (for  conditions 
up  to  2 Km  of  the  atmosphere,  this  is  considered  to  be  valid 

(4) 

assumption) , the  radiative  transfer  equation  can  be  written  as  , 
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A.  /iM  -3.CT)} 
d s t 


where , 


liv)  = the  specific  intensity 
s = integrating  path 
ky  = absorbtion  coefficient 
By(T)  = Planckian  function  at  V and  T 
T = absolute  temperature 

The  solution  of  (1),  evaluated  at  the  ground  is, 

/oo 

I(«)  = +^(s)e^ds 


(1) 


(2) 


where , 


^ = 


ky(s)ds  is  the  optical  depth 


Under  the  assumption  that  thereare  no  sources  lying  outside  of 
the  atmosphere,  equation  (2)  can  be  written  as. 


(3) 


o 

Physically,  equation  (3)  indicates  that  the  received  intensity 

is  the  sum  of  each  incremental  layer  emitting  as  Bj^T)l^s  atten- 

-t 

uated  by  the  intervening  medium  by  C 

Equation  (3)  forms  the  basis  for  all  the  inversion  techniques 
(5) 

that  have  been  used  Mathematically,  equation  (3)  does  not  ad- 

mit to  an  analytical  solution  as  the  variables  involved  cannot  be 
described  analytically.  Consequently,  in  order  to  reach  a meaning- 
ful solution,  some  constraints,  as  determined  by  the  physics  of  the 
atmosphere,  must  be  imposed. 
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B.  Microwave  Radiation  Measurements 


Because  cf  the  simplification  that  can  be  applied  to  equation 
(3),  we  treat  the  microwave  case  first. 

In  the  microwave  region  of  the  spectrum,  the  Planckian  function 
B can  be  simplified  by  the  Ray leigh- Jeans  approximation ( hu  < kT)  as, 

By  ~ ^ ^ /c^  (frequency  units)  (4) 

o 

The  relative  error  in  using  (4)  is  about  0.1%  at  60  GHz  and  300  K. 
Substituting  (4)  into  (3)  we  have  that, 

where , ® 

= brightness  temperature  which  compares  the  received 
radiation  from  a given  direction  to  that  of  a 
blackbody  at  the  temperature 

In  the  troposphere  and  the  microwave  region  considered!  10  to 
70  GHz) , the  principal  atmospheric  gaseous  absorbers  and  hence 
emitters  are  oxygen  and  uncondensed  water  vapor.  The  water  vapor 
absorbtion  is  due  to  a pure  rotational  transition  at  22.2  GHz {1.35cm) 
while  the  oxygen  absorbtion  is  due  to  a band  of  rotational  trans- 
itions staring  at  about  53  GHz  and  extending  to  66  GHz.  An  excell- 
ent article  on  the  spectroscopy  of  the  oxygen  transitions  is  to  be 
found  in  reference  (6)  . 

Before  equation  (5)  can  be  further  simplified,  the  nature  and 
height  dependence  of  the  absorption  coefficient  must  be  known. 
Further,  it  would  be  desirable  to  integrate  equation  (5)  to 
some  finite  limit  rather  than  to  infinity  and  perferably  to  2 KM. 

1.  Oxygen  Microwave  Radiation 

From  the  work  of  reference  (6)  and  (7)  a cut-off  level  can  be 
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established  by  considering  to  what  level  above  which  the  radiation 
emitted  does  not  make  an  appreciable  contribution  to  the  ground. 

This  has  been  done  for  oxygen  in  terms  of  the  ground  absorption 
coefficient  by  reference  (7)  . The  cut-off  level  is  chosen  at  the 
99%  radiation  level;  however,  a more  meaningful  choice  could  be 
made  if  the  exact  experimental  capabilities  of  the  radiometer  in 
use  were  known.  For  example,  if  the  uncertainty  in  the  measurement 
were  .5°K,  then  the  height  level  corresponding  to  this  value  would 
be  chosen.  However,  the  choice  of  99%  will  not  change  the  basic 
conclusions . 

According  to  reference  (7)  , for  a ground  absorption  coefficient 

of  12  DB/KM  or  greater,  the  cut-off  level  is  about  2 KM.  As  the 

(6) 

ground  absorption  coefficient  for  oxygen  is  about  14  DB/KM  , 
we  conclude  that  at  least  to  this  level  we  can  utilize  a constant 
absorption  coefficient  at  the  center  frequency  of  oxygen  band  (60  GHz) 
To  further  simplify  equation  (5),  we  recast  it  as  , 


(6) 


as  <it-  = kds  with  T(s)  expressed  as  T{'t')  and  is  the  cut-off 
optical  depth.  Equation  (6)  can  be  solved  analytically  if  we  make 
the  assumption  that  the  T{'t)  can  be  expressed  as  a Taylor  series 
( this  assumption  would  mean  a nearly  linear  variation  of  the 
temperature  with  V ) . 

Consequently, 

>-  dfT 


TC»)  - T(o^  f ^ ^ ^ ^ • • • 


(7) 


Substituting  into  equation  (6),  it  becomes, 


18 


.ZTc 


-f. 


r^Co)e  dt-  -h 


III  C-fc  ,*01 

^cJ  c»o^ 


/^C 


(8) 


since  kds  and  integrating  the  first  two  terms  we  have, 

TfeCv)  = T>Co)(i  ^ 1 [i  - + o 


within  a given  layer  having  a linear  temperature  gradient,  the 
second  derivative  is  zero  and  (9)  becomes. 


o 


(10) 


If  we  designate  dT/ds  as  a and  choose  to  be  5 or  greater  so  that 

-tc 

C ^ .01,  we  have  that  within  a layer  characterized  by  a linear 
temperature  gradient, 

,11, 


since  we  have  assumed  that  with  a constant  kQ(up  to  2 KM), 

then  tl  = k s defines  the  upper  limit  of  s.  If  we  include  the 
zenith  angle  as  a parameter  then  = k^h^sec  0 . Consequently, 

the  choice  of  0 and  k^  determines  the  height  level, h,  that  the 
measurement  evaluates.  From  equation  (11),  a measurement  of 
that  is  made  for  a ?”c  = 5,  with  the  ground  temperature  T^,  will 
determine  the  slope  a of  the  temperature  gradient.  Assuming  a 
value  of  14  DB/KM  for  the  center  frequency  of  the  oxygen  trans- 
ition and  converting  to  nepers/KM  ( 14/4.3429),  we  obtain  that 
k^  = 3.224.  Consequently,  h^sec  Q = 5/3.224  or  h^sec  9 = 1.55  KM. 
If  the  intial  measurement  is  made  near  the  horizon,  say  5°  elevat- 
ion angle  ( Q =85°) , then  the  effective  radiating  layer  is  about 
155  meters.  It  is  to  be  noted  that  this  is  not  the  spatial  resol- 
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ution  of  the  passive  method  but  the  layer  in  which  a linear  grad- 
ient can  be  ascribed.  For  most  temperature  inversions,  this  is  con- 
sidered adequate  for  accurately  following  the  temperature  gradient. 
However,  if  a finer  resolution  is  required,  then  taking  in  to  acc- 
ount the  earth's  curvature  and  working  a zero  elevation,  spatial 
resolutions  on  the  order  of  15  meters  could  be  achieved. 

To  probe  the  upper  layers,  a second  measurement  is  made  such 
that , 

r^Cv)  .T(o)b-e  t ^L'  Lt,  ^ 

where  'C',  + = "21^  . As  a^^  has  been  determined  from  the  initial 

measurement,  a Z".  can  be  calculated  with  T,  = T + a,h„  of  the 

' 1 o —1  ‘' 

first  measurement.  From  equation  (12) , the  ^is  determined.  Again, 
the  combination  of  k^  and  © is  used  so  that  2]^  = yielding 

the  a in  the  second  layer.  The  process  is  continued  until  the  ent- 
ire 2 KM  is  probed.  From  the  foregoing  analysis  it  is  seen  that  for 
a single  gradient,  only  one  measurement  is  necessary.  For  two  grad- 
ients, two  measurements  are  required  and  so  on.  Based  on  the  fact 
that  10  changes  in  slope  is  about  the  maximum  number  of  changes  that 
occur,  no  more  than  10  measurements  are  necessary  to  characterize 
the  temperature  distribution  of  the  first  2KM  of  the  atmosphere. 

In  general  then,  the  n^*^  measurement  can  be  written  as, 

where  ^ = 5 . 

I 

The  expected  spatial  resolution  will  be  determined  by  the  mea- 
surement errors  and  the  number  of  meaurements  made  but  has  the 
potential,  at  least  in  principal,  of  at  least  15  meters. 
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2.  Water  Vapor  Microwave  Radiation 


The  behavior  of  the  water  vapor  attenuation  in  the  microwave 

region  near  the  resonance  line  (1.35  cm.)  is  remarkable  as  it  can 

(8) 

be  shown  that  the  attenuation  is  independent  of  pressure  and 
depends  only  on  the  fraction  of  water  vapor  present  (i.e.  mixing 
ratio).  However,  even  under  humid  conditions  of  the  atmosphere 
(i.e.  over  ocean  areas),  the  absorption  coefficient  is  not  large 
enough  (->^.66  DB/KM  for  = 21g/M^  at  1.35  cm.)  to  be  apply  the 
method  used  for  oxygen.  Although  the  gross  water  vapor  structure 
could  be  determined,  for  a multi-jagged  variation  where  high  spat- 
ial resolution  is  required,  it  is  not  the  method  of  choice.  Conse- 
quently, it  is  felt  that  the  water  vapor  structure  can  best  be  de- 
termined by  utilizing  the  infrared  water  band  at  6.5^  . 

C.  Infrared  Water  Vapor  Radiation 

In  the  infrared  portion  of  the  spectrum,  the  Rayleigh- Jeans 

approximation  cannot  be  used  and  equation  (3)  must  be  used  as  is. 

Due  to  the  fact  that  the  water  vapor  attenuation  does  not  follow 

(9) 

the  Lambert  law  of  exponential  absorption  , it  is  necessary  to 
cast  equation  (3)  in  a slightly  different  form.  If  the  transmission 
through  a layer  is  t,  then  dt/ds  expresses  both  the  emission  and 
attentuation  from  the  layer.  Consequently, 

dicv)  = BvCs)  a-t  (14) 

For  a layer  of  constant  B(s),  then  the  intensity  received  at  the 
ground  due  to  a layer  at  s is  just, 

Is(^)=  (t.-ti)  (15) 
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where 


= transmission  from  the  bottom  of  the  layer 
t^  = transmission  from  the  top  of  the  layer 
The  measured  intensity  is  just  the  summation  from  each  layer.  Then, 


X(»)  r 6,(0)  (vt.)  +B,(i)Gt,-tC)  • 


(16) 


where, 

B^(0)  = Planckian  function  for  the  ground  temeperature 

t = transmission  at  the  ground  level  equal  to  1 
o 

t^  = transmission  from  the  top  of  the  ground  layer 
t^  = transmission  from  the  top  of  the  second  layer 
B|,(l)  = Planckian  function  for  the  temperature  of  layer  2 
Before  proceeding  to  derive  the  equations  that  determine  the 
spatial  resolution  potential  of  the  infrared  radiation,  an  ex- 
pression relating  the  transmission  to is  necessary.  Although 
the  transmission  cannot  be  expressed  through  an  exponential  att- 


enuation coefficient,  it  can  be  approximated  by 


(7) 


e 


- .b4 


^ Wo/ 


(17) 


where , 

w = optical  thickness  of  water  vapor  in  preciptable  cms. 

w^=  optical  thickness  necessary  to  yield  50%  transmission 

In  addition,  w^  is  a function  of  wavelength  and  can  be  expressed 

as  . 337/  K(X  ^ • The  values  of  K(X  ) have  been  tabulated 

(9) 

by  Green  and  Griggs  . Substituting  the  expression  for  w into 

o 

equation  (17)  we  obtain. 


Now  it  can  be  shown  that  the  equivalent  optical  thickness  of 

water  vapor  in  terms  of  pressure  and  temperature  of  the  atmosphere 

(9) 

can  be  written  as  , _ 

^1. 

Z)..ds 


(19) 


where, 

P = pressure  in  atmospheres 

P = ground  or  local  pressure 
O 

T = absolute  temperature 
T^=  ground  or  local  temperature 
However^^,y can  be  written  as, 


= lO 


M P 


:20) 


RT 

where, 

m = mixing  ratio  (g/Kg) 

R = gas  constant 
M = molecular  weight  of  air 


Substituting  (20)  into  (19),  we  have  that. 


- lO 


w = 


within  an  isothermal  and  isobaric  layer,  this  becomes. 


wts5  - P*Cs)r»'CS>.oS 

■P.V, 


(2i: 


(22) 


Consequently, 


mCs)AS  


Vz 


M Xc"-  P* 

Inserting  the  appropriate  constants,  we  have  that. 


(23) 


23 


(24) 


W(s)aS  - 82. Og?  X »0^  PqT*^  W(AS) 

a«.<)7  P*To''». 


mcs)  AS 


2.  8 32.  X lO* 


P>TH  wCaO 

X ''*- 


As  T (s)  is  known  (oxygen  measurement)  and  P(s)  does  not 
vary  much  and  is  reasonably  predictable  and  with  and  Tq  mea- 
sured quantities,  a measurement  of  w(As)  through  equation  (18) 
determines  the  product  m(s)^s. 

Returning  to  equation  (16),  it  is  necessary  to  limit  the 
intensity  summation  similar  to  way  a cut-off  was  used  in  the 
oxygen  measurement.  The  limit  of  the  summation  process  will  de- 
pend on  the  error  measurement  on  I.  For  the  purposes  of  this  report, 
we  will  assume  that  the  smallest  variation  is  .01  although  the 
actual  value  will  depend  on  the  experimental  capabilities  of 
the  radiometer.  As  B(o)  and  t^  are  known  from  ground  base  mea- 
surements, we  will  choose  a measurement  at  A,  that  produces 
a t^  = .02.  This  in  effect  eliminates  any  term  beyond  B(l)  since 

t ^ t > t and  assures  that  t -t  < .01.  Rewritting  equation  (16) 

1 2 3 2 3 

at  X,  we  have  that, 

I (X.)  = Ct,-tx) 

(25) 

or, 

XU)  = -98 

The  single  measurement  yields  t^.  In  order  to  preclude  the  pos- 
sibilty  that  t^  is  also  equal  to  .02,  we  will  assume  that  the 
mixing  ratio  never  gets  below  40%  humidity  at  any  level  to  2 KM. 
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This  will  correspond  to  a lower  limit  of  about  2 g/Kg.  This  assump- 
tion is  considered  valid  over  the  ocean. 

The  choice  of  t^^  = .02  and  the  ground  measured  value  of  m^, 
produces  a A s at  a given  wavelength  through  equation  (24). With 
the  value  of  t^  obtained  from  equation  (26)  and  the  derived  value 
A s obtained  from  equation  (24),  the  value  of  m^  is  determined. 

The  next  measurement  would  be  made  at  the  Aj^that  produces  a 
t^>.02  (calculated  from  As  and  m^)  and  the  t^  calculated  from  As 
and  m^.  Consequently, 

From  this  measurement,  t^  is  determined  from  which  m^  can  be  cal- 
culated using  the  same  value  of  As  from  t^  and  t2.  This  process 
is  continued  until  the  entire  2KM  of  the  atmosphere  is  probed.  As 
can  be  seen  from  this  process,  the  spatial  resolution  will  be  de- 
termined by  the  number  of  measurements  that  are  made  and  the  limit- 
ations of  the  capabilities  of  the  radiometer. 

In  order  to  assess,  the  potential  spatial  resolution  of  the  pass- 
ive method,  we  compute  the  spatial  resolution  that  could  be  achieved 
from  a humid  ground  layer.  From  equation  (24)  using  = 1 atmosphere 

and  T = 290°K  and  m = 10  g/Kg, 
o 1 

AS  - 9.^12.  X lO^  W Ccm)  (28) 

From  equation  (18),  if  this  layer  is  to  produce  a t^  = .02,  then  for 
K(A)  = 13  (peak  absorption  coefficient  for  water  vapor  occurring  at 
6.6C^)  yields  a w^  = .035  preciptable  cms.  Consequently , AS  = 2874  cm 
or  28.7  meters.  However,  since  As  =Ahsec  & , the  resolution  is  in- 
creased at  large  zenith  angles  so  that  in  principal,  resolutions  on 
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on  the  order  of  2 to  3 meters  ( C7  = 85  ) are  possible.  However, 
as  mentioned  previously,  this  increases  the  number  of  measurements 
that  must  be  made  to  probe  the  2 Km.  of  the  atmosphere.  Additonal- 
ly,  since  the  absorption  coefficient  of  water-vapor  in  the  infra- 
red has  a wide  latitude  (13  to  10~^)  and  in  conjunction  with  the 
zenith  angle,  coarse  or  fine  resolution  measurements  could  be 
achieved. 

Although  the  microwave  radiation  from  atmospheric  oxygen  can 
be  used  to  determine  the  temperature  profile  with  a potential  re- 
solution of  about  15  meters,  it  is  seen  that  the  microwave  radia- 
tion from  atmospheric  water-vapor  does  not  provide  an  attenuation 
coefficient  high  enough  from  which  a high  resolution  profile  could 
be  extracted.  Consequently,  it  is  necessary  to  utilize  the  infra- 
red radiation  from  atmospheric  water-vapor  to  extract  a water- 
vapor  profile. 
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IV.  MICROWAVE  AND  INFRARED  SYSTEM  SENSITIVITY  AND  RESOLUTION 


In  this  section  we  analyze  the  expected  system  sensitivity 
and  resolution.  We  assume  that  the  temperature  profile  has  been 
determined  by  the  microwave  method  and  examine  first  how  many 
meaningful  measurements  can  be  made  using  the  infrared  water- 
vapor  absorption  band  and  what  spatial  resolution  does  this  re- 
present . 

As  we  have  shown  by  the  previous  anlysis,  the  answer  is  com- 
plex and  dependent  on  the  profile  and  the  existing  ground  con- 
ditions. In  addition,  the  spatial  resolution,  as  will  be  shown, 
becomes  non-lnear  which  complicates  arriving  at  a meaningful 
answer.  Although,  using  a spatial  technique  produces  a non- 
linear spatial  resolution,  the  passive  method  is  still  a viable 
method,  as  it  will  be  shown  in  section  VI  that  employing  a grad- 
ient technique  removes  most  of  the  non-linearity.  However,  for 
this  analysis  we  show  how  the  infrared  radiation  from  water-vapor 
will  produce  a non-linear  response  in  terms  of  spatial  resolution. 

A.  Determination  of  Radiometer  Noise 

Before  pursuing  this  analysis,  it  is  necessary  to  consider  the 
limitation  due  to  noise  of  an  infrared  radiometer.  A typical  radi- 
ometer consists  of  a Cassegrain-type  optical  system  which  focusses 
the  radiation  onto  a plane  which  includes  a field-defining  aperture. 
After  mechanical  chopping  of  the  image,  the  radiation  is  relayed  op- 
tically to  an  infrared  detector.  Solid  state  electronics  amplify, 
demodulate  and  filter  the  signal  produced  by  the  detector  producing  an 
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output  voltage  that  is  an  accurate  measure  of  the  received  radiation. 
In  addition,  detectors  are  interchangeable  and  specific  filters  are 
eaisly  interchanged. 

The  radiometric  sensitivity  given  in  terms  of  noise  equivalent 
radiance  (NER)  is  computed  from, 


NEK 


D"  t A 


(1) 


where,  = area  of  the  detector 

^ = output  noise  bandwidth  ( o-  iHz) 

A = area  if  collecting  optics  , 

D*  = detector  figure  of  merit  ( cm (Hz) ^/watts) 
t = transmission  of  the  optical  system 
U>  = solid  angle  field  of  view  (steradians) 

For  a typical  research  radiometer  (e.g,  Barnes  spectral-master.  Models 

12-550  or  12-660) , we  have  that, 

- 3 2 

Ap  = .1mm  X .5mm  ( — 2.5x10  cm  ) 

D*  = 5x10®  (immersed  Bolometer  at  325°K) 

t = .4 

A =90  cm^ 

-4 

CO  = 10  steradians 

M Q O 

Consequently,  for  this  system,  NER  = 10  watts/cm  cO  . We  define 
B(the  blackbody  radiance)  in  terms  of  watts/cm  co  cm^X  where  cm^)i^ 
means  centimeter  wavelength  interval.  If  we  limit  the  radiation 
received  by  the  radiometer  at  particular  height  in  the  atmosphere, 
we  must  have  that, 

B t aX  = nER  (2) 

max 
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where , 


t = transmission  of  radiation  to  the  ground  from  some  level 
of  the  atmosphere 

®max  ~ highest  blackbody  radiance  seen  by  the  radiometer  (i.e. 
highest  temperature) 

= wavelength  bandpass  of  the  radiometer 

If  we  further  restrict  the  radiometer  to  . 1^  wavelength  interval 

at  7yK(i.e.AX=  10  ^cm)  we  have  that, 


t = 10 


-4 


B 


max 


(3) 


For  a normal  lapse  rate,  B is  the  B for  the  ground  temperature. 

O 0 

At  25  C,  B(25°C)  = is  7.16  watts/cm'^CO  at  1 ju  . Consequently, 

“ 5 

t = 10  . This  indicates  that  the  radiometer  sensitivity  is  on  the 

-5 

order  of  10  B . However,  this  is  necessarily  too  sensitive  and  we 
max 

assume  that  by  a suitable  change  of  one  of  the  parameters  in  equation 

-4 

(1),  we  can  degrade  the  sensitivity  to  10  B . We  will  call  this  the 

max 

noise  equivalent  blackbody  radiance  (NEB) . Consequently,  any  signal 
received  can  be  written  as. 


-4 

Signal  = SNR  10  B 


max 


(4) 


where, SNR  is  the  signal  to  noise  ratio. 


B.  Determination  of  Sensitivity 

According  to  section  III, part  C,  the  measured  intensity  or  effective 
B is  given  by, 

I/ZS^K  = B^(l-t^)  + B2(t^-t2)  + •••  Bn^^n"^n+1^ 

We  rewrite  equation  (5)  as. 
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1 


I/aX  - = B2(t3^-t2)  + 

where,  summation  of  all  the  upper  layers.  If  t^ 

that  the  ^i.  *=  0 (i.e.  t^  = 10  then, 

2 


I/aX  - Bj^d-tj^)  = B2(t3^-t2) 


The  magnitude  of  the  signal  is  given  by  equation  (4) . 


= SNR  10  '^B 


max 


or  dividing  by  t^  we  have  that. 


^2  ®2 

SNR  = _f { t^/t^  - 1) 


-4 

10  B 


max 


-4 

since  t^  = 10  , 


SNR  = 


( t^/t^  -1) 


max 


i.e.,  independent  of  t2- 

Before  proceeding,  we  require  a known  temperature 

first  case,  we  choose  a linear  gradient  starting  with 

o 

ground  and  decreasing  to  288  K at  2 KM.  For  this  case 
SNR  becomes , 

B(288) 


B(298) 


or, 


<S™>min  = •”<tlA2  - 1> 


From  section  III, part  C,  we  showed  that. 


rv.55  .55 
-1 . 90Ksec  y w 


= e 


where  0 = the  zenith  angle.  Hence, 


(6) 

is  chosen  such 

(7) 

Consequently , 

(8) 

(9) 

dO) 

profile.  As  a 
298°K  at  the 
, the  minimum 

(11) 

(12) 

(13) 
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n.55  .55  .55 

, -1.90Ksec  0 (w,  - w.,  ) 

t /t  = e 12 


(14) 


We  will  call  Ksec 


e-”  , , 


, -1,90  6 (w. 

t^/t^  = e 1 


.55 


and  write  that, 
.55 

- W2  ) 


(15) 


In  order  for  equation  (7)  to  be  valid,  we  also  must  have  that, 

55 


t2  = 10  = e 


-4  -1.90  ^ w^ 


(16) 


or. 


$ = in 


(10"‘) 


1.90  w. 


. 55 


(17: 


Consequently,  a given  W2  defines  ^ . This  indicates  that  any 
combination  of  K and  sec  0 produces  a given  ^ . For  the 

purposes  of  this  analysis  we  will  assume  that  sec  0 varies  from 

1 to  7 corresponding  to  0 = 0°  to  0 = 82°.  For  © = 82°,  we  can 

still  assume  a flat  earth.  The  variation  of  K aoes  from  a peak  of 

about  13  at  6.70  yu  and  we  chose  a lower  limit  of  K=1  at  about 

7.50^  . Consequently  ^ can  vary  from  1 to  about  38. 

A plot  of  versus  W2  is  shown  in  figure  3 according  to  equation 
(17).  On  the  same  curve  is  plotted  the  expected  SNR  for  a given 
( ^ . This  was  computed  by  combining  (12)and  (15). 

These  curves  indicate  the  expected  SNR  for  a given  difference  of  w 
in  a two  layer  system. 

It  is  to  be  noted  that  in  this  method,  we  are  always  working  in  a 
wavelength  region  where  the  radiation  is  optically  black.  Consequently, 
the  measurements  in  intensity  at  different  , will  only  reflect  the 
changing  B or  the  ef fective^B/t 

Before  proceeding,  we  now  require  the  ground  value  of  the  water  vapor 
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densityy^n/  or  its  equivalent  the  mixing  ratio.  In  addition,  we 
choose  the  spatial  resolution  required.  The  combination  of  measured 
^yOivand  aH  produces  , the  optical  depth  of  at  A H .Con- 

sequently, for  a given  $ , the  term  (1-t^)  in  equation  (7)  is 
defined  which  in  turn  defines  B2  (ti-t2) . 

The  determination  of  $ to  be  used  to  calculate  B2  (ti~t2)  in 
equation  (7)  is  obtained  as  follows: 

The  or  w at  any  given  increment  AH  must  lie  between  Wg 

( w at  saturation)  and  some  limiting  value,  say  20%  humidity  or  .2  w . 
This  is  indicated  in  figure  4.  Plotted  between  them  is  the  w for  the 

linear  temperature  profile  given  above  and  an  exponential  water  vapor 

3 3 

profile  starting  at  20  g/M  at  the  ground  and  ending  at  10  g/M  at  2 

KM.  This  w curve  was  calculated  at  50  meter  increments.  We  use  this 

curve  as  a model  to  the  determine  the  spatial  resolution. 

At  a A H = 50  meters  and  a the  w at  50  meters  is  .10 

pr.  cm.  If  we  continue  at  50  meter  increments,  the  value  of  W2  at  100 

meters  for  a saturated  layer  is  .3  pr.  cm.  while  for  a 20%  humidity 

it  is  .12  pr.  cm..  Hence  expect  that  the  measured  value  between  50  and 

100  meters  should  fall  between  the  values  computed  as  follows: 

Wg(50  to  100)  = Wg(lOO)  - Wg(50)  = (.3  -.15)  = .15 
W2g(50  to  100)=  W2q(100)  - W2q(50)  = (.06  -.03)  = .03 

Consequently  for  the  50  to  the  100  meter  layer, 

Wp  = .l(+  .15]  .25 

Referring  to  the  curve  of  figure  3,  this  produces  a range  of 
^ = 15.0  to  ^ = 10.2.  The  measured  value  of  I/^X  - Bj^d  -tj^)  can  be 
obtained  as  is  scanned  from  15.0  to  10.2.  Each  value  of  ^ will  pro- 
duce a signal  at  a given  SNR.  However,  the-S  values  for  saturation  and 
20%  humidity  are  .184  and  .043  respectively.  At  - S = .184  at  10.2, 
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the  SNR  is  27  while  at  -6=  .043  and  $ = 15.0,  the  SNR  is  2.  Since 
the  actual  value  of  w is  .20,  this  would  produce  an  SNR  of  about  14 
at  ^ = 11.8.  Scanning  the  ^ from  15.0  to  10.2  would  produce  the  SNR 
that  is  consistent  with  the  actual  W2. 

Once  VI 2 is  determined,  the  process  is  continued  for  the  third  layer. 
Now  tj^  represents  the  transmi  ssion  to  the  ground  from  the  top  of  the 
first  two  layers  and  the  t2  from  the  top  of  the  third  layer.  In  each 
layer,  the  limits  of  saturation  and  20%  humidity  determines  the  range 
of  ^ that  is  to  be  scanned  and  the  SNR  obtained  that  is  consistent  with 
the  - S for  each  ^ results  in  the  value  of  w for  that  layer. 

To  obtain  a measure  of  the  expected  resolution  for  this  process 
and  the  number  of  measurements  that  can  be  made,  we  examine  the  process 
in  reverse  and  calculate  the  minimum  thickness  of  the  atmosphere  that  is 
necessary  that  will  produce  a SNR  = 1,  starting  from  the  top  of  the 
2 KM  layer  using  the  calculated  w curve  for  the  model yO  profile. This 
is  computed  from, 

^ _ ln(2.27) 

OsNR=l 

-1.90 


and  from. 


w 


= ( S+  ^2 


55 


1.818 


(19) 


The  results  are  shown  in  Table  I.  From  this  we  conclude  that  for  a 
nearly  satured profile,  neglecting  the  first  100  meters,  we  can 
make  about  17  measurements.  However,  the  spatial  resolution  is  non- 
linear, being  fine  at  the  bottom  of  the  atmosphere  and  becoming 
coarser  as  higher  layers  are  probed.  As  the  course  of  radio  waves  are 
affected  by  w rather  than  , we  believe  that  the  coarse  resolution 
at  the  too  layers  is  not  a serious  restriction  of  the  method. 
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TABLE  I* 


H (Meters) 

'"2 

$ 

S 

w -55 
'^1 

. 55 

Wi 

2000 

2.96 

2.70 

-.160 

1.66 

2.51 

1600 

2.51 

2.92 

-.148 

1.51 

2.12 

1300 

2.12 

3.20 

-.135 

1.38 

1.79 

1100 

1.79 

3.50 

-.123 

1.25 

1.51 

850 

1.51 

3.80 

-.114 

1.14 

1.27 

700 

1.27 

4.25 

-.102 

1.04 

1.07 

580 

1.07 

4 . 70 

-.092 

. 946 

.904 

480 

.904 

5.10 

-.085 

.861 

.762 

400 

.762 

5.60 

-.077 

.784 

.643 

330 

. 643 

6.20 

-.0696 

.715 

. 543 

280 

. 543 

6.70 

-.0644 

.650 

.457 

230 

.457 

7.40 

-.0583 

. 592 

. 385 

200 

. 385 

8.20 

-.0526 

. 539 

. 325 

170 

.325 

9.00 

-.0479 

. 491 

.274 

140 

.274 

9.90 

-.0436 

.447 

.231 

120 

100 

1 - 1 

.231 
. 195 

10.80 

-.0400 

.407 

. 195 

*Tabulation  of  W values  for  the  mode]  profile  used  to  determine  the 
expected  resolution  (see  text  above) . 

C.  Microwave  System 

The  microwave  radiometer  spatial  resolution  and  temperature 
sensitivity  is  a function  of  the  requisite  parameters  for  a via- 
ble probe.  For  example,  as  the  normal  temperature  gradient  of  the 
atmosphere  is  about  -6.5  °K/Km  (-.0065  °K/Km)  and  a spatial  reso- 
lution of  50  meters  is  required,  then  the  microwave  radiometer 
should  be  able  to  resolve  about  .3°K. 

The  temperature  sensitivity  of  the  radiometer  can  be  expressed 


as , 


where. 


CTe 

° /hi 


(20) 


T . = root-mean  square  tonperature  fluctuation (minimum  detectable 

tenperature) 

N = noise  figure 

t = integration  time 

Tq  = ambient  temperature 

C = constant  between  1.5  and  2 (depends  on  mode  of  operation) 

T^  = excess  receiver  noise  tenperature 


For  a sytem  operating  at  60  C3]z  (oxygen  band) , with  an  IF  bandpass  of  20  MHz 


(typical) , the  T . is  a function  of  the  excess  receiver  noise.  Assuming  an  in- 
mm 

tegration  time  of  10  sec.  and  a T^  of  290°K  with  C = 2,  we  have  that. 
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System  Noise  (DB) 

Sensitivity  (°K) 

20 

4 

17 

2 

10 

.37 

8 

.22 

5 

.09 

Consequently,  a system  noise  of  about  10  DB  is  necessary  to  meet 
this  requirement. 

In  terms  of  the  spatial  resolution  of  the  radiometer,  the  gain 
of  the  antenna  used  determines  the  area  over  which  the  received 
energy  is  integrated.  Hence  the  beamwidth  is  a measure  of  the  re- 
solving power  of  the  antenna.  For  a paraboloid  antenna,  the  eff- 
ective area.  A,  is  about  55%  of  the  physical  area.  A general  re- 
lation of  the  antenna  gain  G,  is  given  by, 

Ct  - A (21) 

X" 

At  60  GHz  and  a diameter  of  the  antenna  of  1 foot,  the  gain 
(55%  efficiency)  is  about  40  DB.  The  angular  resolution  of  the 
antenna  is  given  by. 


where  X = wavelength  of  the  received  radiation (. 5cm) 

D = effective  diameter  of  the  antenna  (D^  = Dv/.  55  ) . 
o o 

Consequently,  the  angular  resolution  is  about  .022  radians ( ~1 . 3°) . 
At  2 KM  (limit  of  the  temperature  probe),  this  represents  a 
physical  area  of  about  44  meters  in  diameter  which  is  consistent 
with  the  assumed  temperature  spatial  resolution  requirement. 

In  practice,  the  antenna  gain  is  determined  empirically.  For 
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the  stated  requirements,  a 1 foot  diameter  antenna  with  a noise 
figure  of  10  DB,  appears  to  meet  the  needs  of  a viable  probe. 

Both  system  cost  and  complexity  cannot  be  determined  but  must 
be  evaluated  when  the  results  of  the  error  propagation  analysis 
is  made. 

As  microwave  radiometry  is  a well  developed  technology,  it  is 
concluded  that  a properly  designed  system  can  adequately  meet  any 
necessary  requirements. 


D-  Infrared  System 

The  spatial  resolution  of  an  infrared  system  utilizing  a 
Cassegranian  telescope  as  the  receiver  more  than  adequately 
meets  the  spatial  resolution  requirements  (due  to  much  shorter 
wavelength  ~ 6 . 5^ ) . The  sensitivity  of  the  system  to  water  vapor 
concentration  changes  will  be  a function  of  the  limits  imposed 
by  the  error  analysis,  the  requisite  requirements  of  the  probe 
and  the  ability  of  the  radiometer  to  measure  intensity  changes. 

The  noise  equivalent  radiance  (NER)  of  the  radiometer  can  be 
computed  from. 


NER 


X)"e  A CO 


(23) 


where,  CX,  = detector  area 

Af  = output  noise  bandwidth 
A = area  of  the  collecting  aperture 
OJ>  = solid  angle  field  of  view 
D*  = detector  figure  of  merit 
e = transmission  of  the  optical  system 

For  a typical  system,  we  have  that  CX  = .12mm  x .48mm,  Af  = 1 Hz, 

2 8 ^ 

A = 8.5  cm  , to  = .0012  sr,  D*  = 10  cm  ([}z ) ^/watt  and  e = .4.  Using 

“8  y 

these  values  NER  = 6x  10  watts/cm  -sr. 

To  translate  this  into  an  equivalent  sensitivity  for  water 
vapor  concentration,  we  proceed  as  follows; 
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For  a given  radiating  layer,  the  emitted  intensity  is  I = 

1j 

Bj^(t2  - t^)  where  the  symbols  were  defined  in  progress  report  #2. 

As  t2  is  the  transmission  from  the  top  of  the  layer  to  the  ground, 
it  can  be  expressed  as  t2  = where  tj^  is  the  transmission 

through  the  layer.  Then  = (tj^  - tLtj^)Bj^.  A change  in  the  emitted 
intensity  due  to  a change  in  the  transmission  of  the  layer  from  a 
change  in  the  water  vapor  concentration  and  assuming  a constant 
temperature  in  the  layer  results  in 

shown  in  section  III,  part  C chat  tj^  is  related  to  the  optical 
thickness  w through. 


t e 


— 1.90  C>^)  w 


(24) 


a change  in  t^  due  to  a change  in  w becomes. 


.55 


^b,  ^ -190KCX)e  ( 5-5V  (25) 


In  terms  of  the  mixing  ratio  m,  the  w is  in  turn  related  to  m 
through  mAs  S'  8.50x  lO^w,  assuming  a = 1 atmosphere  and  a 
= 300°K.  Using  a K ( >. ) = 13  (peak  of  the  water  vapor  absorbtion) 
and  substituting  into  equation  (25)  we  obtain  that, 


„ _ 8.50  X I O-^  lu 


(26) 


If  we  assume  a spatial  resolution  (in  depth)  of  50  meters  as 
before  and  a typical  value  of  w = .01  pr . cm.,  with  t^  nearly  1, 


we  have  that. 


A no  - 


O A Tl 


(27) 


We  equate  the  A I to  the  sensitivity  of  the  radiometer  and 

Li 

obtain  that  the  sensitivity  in  terms  of  the  mixing  ratio  is, 


Arr.  = 


(28) 
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The  is  the  blackbody  radiation  at  the  bandpass  and  wavelength 
of  the  radiometer.  At  6.5  microns  with  T = 300^K  and  assuming  a 
bandpass  of  0.1  microns  for  the  radiometer,  = 4xl0~^  watts/cm^- 
steradian.  Substituting  the  calculated  value  of  the  NER,  we  ob- 
tain that, 


9/K3 


(29) 


Although  this  is  an  approximate  value,  it  indicates  that,  poten- 
that  the  method  has  a sensitivity  that  exceeds  any  reaso- 
nable requirements  for  a probe. 
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FIGURE  3.  CURVES  USED  TO  DETERMINE  THE  SNR  FOR  A GIVEN  W , d)  AND  <J  • (SEE  TEXT) 


(Meters) 


FIGURE  4 . CALCULATED  VALUES  OF  W FOR  THREE yO  PROFILES  VERSUS  H 

(SEE  TEXT) 
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V,  EFFECT  OF  VISIBILITY  ON  PASSIVE  PROBES 


In  order  for  a probe  for  water-vapor  and  temperature  to  be 
useful,  whether  active  or  passive,  it  must  be  capable  of  opera- 
ting under  varying  degrees  of  visibility.  The  effect  of  visual 
range  (i.e.  visibility)  of  1 to  10  Km  on  both  the  infrared  and 
microwave  passive  methods  is  evaluated  in  this  section. 


As  visual  range  is  a subjective  parameter  while  the  proposed 
passive  methods  are  objective  measurements,  a criterion  must  be 
established  to  relate  the  two.  Visual  range,  as  defined,  depends 
on  the  ability  of  an  observer  (i.e.  the  "standard"  eye)  to  dis- 
tinguish an  object  from  it's  bac)<:ground  light  (contrast).  In  its 

most  general  form,  the  reduction  in  contrast  for  an  observer,  un- 

(11) 

der  different  s)cy  brightness  conditions  is  given  by  , 

(1) 

where,  B^  and  B^  are  the  inherent  luminances  of  the  object  and 

bac)cground  respectively,  B and  B'  are  the  reduction  in  luminance 

K R 

at  the  distance  R of  the  object  and  bac)<ground  respectively,  <3^,  the 
extinction  coefficient  for  a homogeneous  layer  of  the  atmosphere 
and  R,  the  optical  range;  defined  as  the  equivalent  horizontal  dis- 
tance which  produces  the  same  attentuation  as  the  actual  range  R. 
Calling  the  inherent  contrast  C^,  we  have  that. 
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(3) 


and  the  apparent  contrast  Cj^  as, 


r - ^ 

' — ^ ^ 

j5r 


we  have  that  equation  (1)  can  be  written  as, 

Cr  = Co  (Bo/E>p.)e  ° (4) 

Due  to  the  fact  that  the  eye  contrast  threshold  ^ , depends  on 
the  angular  subtense  of  the  object  and  that  for  a particular  ob- 
ject this  in  turn  depends  on  the  range  and  that  the  resulting 
contrast  depends  on  the  luminance  of  the  background,  equation  (4) 
admits  to  a range  of  values.  Nomograms  for  a variety  of  luminance 

values  and  object  sizes  are  available  to  determine  visual  range 

( 12 ) 

for  a given  set  of  conditions  . However,  for  the  purposes  of 

evaluating  the  passive  methods,  we  utilize  equation  (4)  for  the 

case  of  an  object  seen  against  the  horizon  sky.  If  C is  greater 

o 

than  , then  equation  (4)  becomes. 


€ - 


Ce 


-C3;v 


(5) 


where,  V = visual  range 

luminal  contrast. 

As  we  are  concerned  with  the  evaluation  of  C5^ , we  assume  a black 
object  {C^  = -1)  and  equation  (5)  can  be  rearranged  to. 

Before  CSgCan  be  evaluated,  a value  for  €:  must  be  assigned.  If  we 
use  the  meteorological  range  as  the  criterion  for  visibility,  then 
^ = .02  and  equation  (6)  can  be  written  as. 
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(7) 


Oo  ^ 3.912.  A/;^ 


where,  V2  = the  meteorological  range.  Consequently,  a meteoro- 
logical range  of  1 to  10  KM  produces  a value  of  CToOf  3,912  KM  ^ 
to  .3912  KM  This  range  of  extinction  coefficients  must  now  be 
related  to  some  power  law  with  wavelength  in  order  to  assess  its' 
effect  in  the  infrared  and  the  microwave  regions. 

The  most  general  form  of  CT  can  be  written  as, 

V o \ 

CT  = C,  A t A (3) 


where,  c^^  and  are  constants,  X the  wavelength  and  the  second 

term  reflects  the  Rayleigh  extinction  coefficient  for  pure  air. 

(11) 

As  for  pure  dry  air,  the  visual  range  is  over  350  KM  , we 
are  justified  in  considering  only  the  first  term  in  equation  (8) 
when  dealing  with  a visual  range  of  1 to  10  KM.  Consequently, 

O’  = C,  ^ where  the  exponent  n must  be  evaluated  on  the  basis 
of  the  type  of  scattering  and  absorbing  particle  involved. 

In  order  to  evaluate  the  exponent  n,  recourse  is  made  to  the 
Mie  theory  of  scattering  in  which  the  cross-section  for  extinc- 
tion can  be  wriiten  as. 


O'  = Ni  Trr"-K 


(9) 


where,  N = number  of  particles  cm~^ 
r = radius  of  the  particle 
K = efficiency  for  extinction 

K is  ratio  of  the  extinction  coefficient  to  the  geometric  cross- 
section  of  the  particle.  From  the  Mie  theory,  the  K is  a function 
of  <X,  the  size  parameter,  equal  to  2Trr/X.  For  a dielectric  sphere 
(i.e.  a water  drop)  with  a refractive  index  of  1.33,  K is  about  4 
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when  OC  = l.  Above  OC  = 1,  the  K fluctuates  about  K = 2 and 
approaches  2 for  large  values  of  CX  and  eventually  approaches  1 
in  the  geometric  region.  From  equation  (8)  and  equation  (9)  we 
can  write, 

c.^r/x) 

(10) 

* 

Taking  the  log  of  equation  (10)  yields, 

1 oc^  K +-  I o<^(N  it  Ir^)  - I ocj  C|  - n 1 ocj  ( 

Consequently,  for  a given  value  of  r and  X , K can  be  calculated 
from  Mie  theory  for  a given  refractive  index,  allowing  the  deter-  | 

mination  of  n.  This  has  been  done  for  water  spheres  (m  = 1.33)  and  | 

completely  absorbing  spheres  (m  = c>o  ) by  reference  (11).  For  water  j 

spheres,  the  value  of  n = -4  for  very  small  spheres  (Rayleigh  ex-  | 

tinction) , becomes  n = -2  for  r about  ^2  and  n = 0 for  r = X . ■ 

It  then  becomes  temporarily  positive  and  oscillates  about  zero  as 
r/X  increases.  For  completely  absorbing  spheres,  the  transition 
from  n = -4  to  n = 0 occurs  at  r/X  = .2.  Although  the  Mie  theory 
predicts  an  oscillating  value  for  n for  large  values  of  ^ , these 
effects  are  not  observed  as  atmospheric  aerosols  are  mainly  poly- 
disperse  resulting  in  a dampning  of  the  oscillation  in  n.  Consequen- 
tly, depending  on  the  value  of  m,  a critical  value  of  r is  reached 
where  the  extinction  becomes  independent  of  wavelength. 

Accordingly,  the  extinction  coefficient  for  atmospheric  aerosols 

can  be  used  to  determine  the  value  of  n.  Empirically  it  has  been 
(13) 

found  that  , 


Although  this  equation  has  been  evaluated  for  a range  of  1 to  5 KM, 
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for  the  purposes  of  this  calculation,  we  will  assume  that  it  is 
also  valid  up  to  10  KM.  Combining  equation  (7)  with  (9)  we  have 


that , 


<N> 


3. 912 


(13) 


where  the  brackets  indicate  the  average  value.  As  equation  (12) 
permits  a determination  of  n which  in  turn  defines  both  K ( oC ) and 
r,  ( using  and  average  value  of  X = . ^ for  the  visible  ),  permits 
and  estimation  of  <^N)>from  equation  (13). 

Assuming  that  water  droplets  produce  the  reduction  in  visibility, 
(reasonable  assumption  for  open  ocean  areas) , a range  of  1 KM  pro- 
duces  an  <,r'>  = . 43^  and  <^n'>  2.2x  10  /m^,  corresponding  to  an  n = 

.59  and  K "=  3 . For  a visual  range  of  10  KM,  n = 1.26,  r/X«.65,<r>  = 
.33^,  K = 3 and  <Cn'>=  3.8x  10  /m^.  Consequently,  we  are  to  assess 

Q -3  9 3 

the  effects  of  water  droplets  in  the  range  of  4x  10  /m  to  2x  10  /m  , 
varying  in  size  from  .4  to  .3ywfor  the  infrared  and  the  microwave 
regions . 

Before  addressing  this  problem,  it  is  instructive  to  compare 

the  obscuring  power  of  rainfall  with  that  of  haze  considered  above. 

(11) 

It  can  be  shown  that  , 

(5  -5.2X  IO"^X^(cm-9 

where  X = number  of  raindrops/cm^-sec  and  the  drops  are 

enough  for  K to  be  2 (greater  than  10 X).  From  equation 

noting  that  Z,  the  precipitation  rate  (cm/sec.),  is  XVj. 

is  the  volume  of  a single  raindrop,  we  have  that, 

Z ( ^^/sec)  S .75  X (O^- 

V,  Ccn) 

Assuming  a raindrop  size  of  1 mm  radius,  we  have  that. 


(14) 

large 
(7)  and 
where  v^. 

(15) 
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(16) 


ZC^'^/sec)  ^ 


Ci:;  -TT  A \0' 


VjCcm) 


5 -2 

Consequently,  for  a of  1 KM  (10  cm),  Z = TT  x 10  and  for  V2  = 

10  KM,  Z =Tr X 10”^  cm/sec.  As  rainfall  seldom  exceeds  2 to  3"/hr., 
we  conclude  that  rainfall  can  have  less  obscuring  power  than  haze 
(i.e.  greater  visual  range).  These  results  will  be  used  when  consid- 
ering the  infrared  and  microwave  regions. 

In  order  to  assess  the  effect  of  haze  in  the  infrared  and  micro- 
wave  regions,  we  utilize  the  general  form  of  the  extinction  co- 
efficient from  the  Mie  theory.  For  an  absorbing  particle  with  a 
complex  refractive  index  ( water  in  the  infrared  and  microwave) , 
the  scattering  coefficient  for  values  of  CX.  < 1 is  given  by. 


3Tr 


(17) 


where  ]t  =(m^  - iVon  + 2).  The  absorbtion  coefficient  is  given  by. 


Qa 


Cs^  ^ I (-  k} 


TT 


(18) 


7 2 

where  Im  is  the  imaginary  part  of  -(m^  - l)/(m  + 2) . Equations 

(17)  and  (18)  combined  produce  the  extinction  coefficient  per 
particle  and  can  be  related  toCTof  the  visual  range  by  multiplying 
by  calculated  above. 

For  water  droplets  in  the  infrared,  the  complex  ref reactive 
index  at  6^ is, 

nn-  a 1,304  - .1  i 


For  water  spheres  with  <r>  = .35^,o<.^  .37.  From  equation  (19), 
(m*)^  = 1.69  - .26i  and  k = .18  - .08i  with  = .04.  Accordingly, 
Qg  = 7 . 8x  10  m^  and  = 4.64x  10~^^  m^ . Consequently,  Q =4.7 
X 10  m . From  the  visual  range,  the  maximum  was  found  to  be 
about  2x  10^/m^.  Hence  C = 9.4x  10  m ^ or  CS'  = .094  KM~^.  In 
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terms  of  optical  depth  used  in  section  III,  part  A,  even  at  2 KM, 
the  effect  on  the  ^ due  to  water  droplets  is  not  significant  com- 
pared to  the  Z^of  water  vapor.  Consequently,  we  conclude  that  even 
with  a visual  range  of  1 KM,  the  effect  on  the  infrared  measurement 
is  not  significant.  In  the  case  of  rainfall,  the  effect  would  be 
about  the  same  as  in  the  visible  asa>10  and  the  extinction  coeff- 
icient is  nearly  the  same  for  absorbing  and  non-absorbing  particle. 

For  the  microwave  region  for  haze  particles,  as  the  value  of 
-4 

is  so  small  { ~ 10  ) , only  Rayleigh  scattering  prevails  and  the 

effect  on  (3”  of  the  oxygen  band  (^5mm)  is  even  less  than  the  infra- 
red region.  However,  for  rainfall,  where  the  drop  size  is  comparable 
to  the  microwave  wavelength,  equations  (18)  and  (19)  can  be  used. 

The  values  for  |ld  and  Im  (-k)  for  water  at  .62cm  has  been  tabula- 

ted^^^^and  the  attentuation  coefficient  per  mm  rainfall/hr.  has 

(15) 

been  calculated  by  Medhurst  . At  60  GHz  ( oxygen  frequency) , the 
maximum  rainfall  attenuation  coefficient  measured  is  about  .85  DB/KM 
per  mm/hr.  From  the  previous  calculation  utilizing  the  visual  range 
it  was  estimated  that  a visual  range  of  10  KM  corresponds  to  a rain- 
fall of  about  10  cm/hr., (a  rate  only  occuring  for  a short  time  at 
frontal  regions).  If  we  assume  a rainfall  of  1 cm/hr.  (correspond- 
ing to  a heavy  downfall) , then  it  is  expected  that  there  will  be  an 
appreciable  effect  on  the  microwave  method  as  the  ground  value  for 
oxygen  absorbtion  is  about  14  DB/KM. 

It  is  concluded  from  the  above  analysis  that  for  visual  ranges 
of  1 to  10  KM  due  to  haze  or  fog,  the  infrared  and  microwave  methods 
are  not  significantly  affected  and  will  not  result  in  appreciable 
errors;  however,  for  rainfall,  the  microwave  region  is  affected  more 
than  the  infrared  and  corrections  are  necessary  to  utilize  the  method. 
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VI.  RECOVERY  OF  TEMPERATURE  AND  WATER-VAPOR  PROFILES  FROM  THE 

PASSIVE  PROBES 


Rather  than  examining  every  combination  of  temperature  and 
water-vapor  distributions,  which  would  detract  from  the  main 
goals  of  this  study,  we  examine  in  this  section,  a temperature 
distribution  with  two  inversions  combined  with  a water-vapor  dis- 
tribution with  five  inversions  up  to  two  kilometers.  The  rationale 
for  this  combination  stems  from  the  belief  that  this  combination 
represents  a worse  case  and  if  these  distributions  can  be  reco- 
vered satisfactorily,  then  any  combination  having  less  inversions 
will  also  be  recovered  satisfactorily. 

The  temperature  profile  used  is  depicted  in  figure  5 and  the 
water-vapor  profile  used  in  combination  is  shown  in  figure  6. 
Although  this  may  not  represent  a real  case,  the  results  should 
be  an  indication  of  the  enormous  potential  of  the  passive  method. 


A.  Temperature  Profile  Analysis 

As  mentioned  in  the  previous  sections,  the  measured  temperature 
by  a microwave  radiometer  at  the  ground  due  to  a given  distribution 
is  given  by. 


J 


(1) 


where,  is  the  ground  temperature,  A.h  is  the  increment  of  height, 
<Pois  k^sec  0 , k^is  the  ground  value  oxygen  attenuation  coeffi- 
cient, © is  the  zenith  angle  and  T^  is  the  temperature  of  each  suc- 


ceeding layer.  If  we  let  x = e 


then  we  have  that. 


49 


PRECSDINQ  PADS  BUNK-NOT  flU'lED 


(2) 


7;,=  6-x)  I^To  (To+-aAb)e'^‘>''V(T«  + 2aAK)e^'^^^^^ 


where,  - t + naAh  ( linear  gradient  within  each  layer).  Rearran- 


ging (2)  we  have  that. 


Tv,  = 6-x)j^To  4- (To qa1-))x  +•  (1^+ 2aAh)xV  • • -J 


(3) 


or. 


T„  =(i-x) 


To  I I + X +x’*’+x’+  - j + QaK^X-H  2x*-4  5x’4  - 


(4) 


or. 


T„  - (l-x)  [to  { I 4X-fx'4.  . x'^}4-  aXAh^l  42X  43x^4... NX^'H 


(5) 


If  we  let  N go  to  infinity,  we  have  that. 


T„  = (1  -a)  |to  ^ 


or. 


TTi  = To  4-  ax.ah 
0-x) 

Substituting  back  for  x,  we  have  that, 

- 4‘o'^b  I 

-r-T-i-®-® 

' M ” 'O 


(1  - 


(6) 


(7) 


(8) 


Using  L’Hospital’s  rule  and  letting^hgo  to  0,  we  finally  obtain  that. 
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Tn  = '^o  * 


o 


(9) 


This  represents  the  integrated  value  of  measured  temperature.  It  is 
interesting  to  note  , that  this  same  result  was  obtained  in 
section  Ill/squation  (11)  using  an  entirely  different  approach. 

If  the  remains  finite,  then  the  equivalent  to  equation  (9)  is 
given  by, 

+ Q.>hx 

(.  I-X 


where  N represents  the  number  of  layers  taken  with  a given  ^h. 

We  have  programmed  both  equations  (9)  and  (10)  for  the  temperature 
profile  given  in  figure  1.  We  have  compared  the  for  ^h=o(true 
measured  temperature)  with  A h =50 , 100  and  200  meters.  The  results  of 
these  computations  are  not  reproduced  in  this  report. 

The  integrated  case  results  ( i.e. what  is  measured  by  the  radio- 
meter) is  shown  in  figure  7 for  the  temperature  profile  of  figure  5. 

The  actual  values  of  ^^used  for  the  computation  were: 

(1)  .030/meter 

(2)  .021/meter 

(3)  .010/meter 

(4)  .005/meter 

(5)  .003/meter 

(6)  .001/meter 

For  the  .001/meter  case,  it  was  necessary  to  integrate  up  to  6.4  km. 
using  a constant  gradient, 2 km  on.  This  point  was  excluded  from  figure?. 
For  the  .003  case,  the  integration  went  to  3.6  km  before  reaching  a 
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constant  value.  For  the  .005  case,  it  stabilized  at  about  2 km  with  the 
other  cases  stabilizing  at  less  than  2 km. 

From  the  results  of  these  computer  calculations,  we  were  able  to 
determine  the  cut-off  height  over  which  a given  is  effective.  This 
is  shown  in  figure  8.  It  is  to  be  noted  that  for  values  of  less 
than  .0025,  the  cut-off  height  is  greater  than  2 km.  Consequently,  in 
practice,  we  would  utlize  values  of  equal  to. 002  or  greater. 

According  to  equation  (9)  , a measurement  of  at  given  will 
result  in  an  apparent  a.  We  have  calculated  the  apparent  ^ as  a 
function  This  is  shown  in  figure  9.  As  noted  from  the  figure,  the 

apparent  ^ remains  constant  until  about  4^ = .007.  From  figure  8 the 
effective  cut-off  is  700  meters.  Consequently,  we  can  conclude  that 
above  700  meters  there  is  a change  in  gradient. 

In  order  to  recover  the  profile  of  figure  5,  we  would  conclude  that 
up  to  700  meters,  the  gradient  is  constant  and  changes  above  700  meters. 
The  gradient  for  the  first  700  meters  is  computed  from  equation  (9). 
Consequently,  any  temperature  measurement  from  4)^ equal  to  .021  to 
.007  would  yield  a gradient  of  - . 00 5°K/meter . In  order  to  recover  the 
rest  of  the  profile  we  note  that,  for  a temperature  distribution 
with  inversions  we  have  that, 


(11) 


Since  the  apparent  gradient  is  given  by. 


(12) 
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FIGURE  8.  EMPIRICALLY  DETERMINED VERSUS  (^EE  TEXT) 
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Apparent  Gradient  a (°K/ineter) 

FIGURE  9.  APPARENT  GRADIENT  CALCULATED  FROM  MEASURED  TEMPERATURE 
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we  have  that 


a - Q I 


o,  e 


Q^e 


(13) 


where  the  represent  the  heights  at  which  the  inversions  occur.  To 
determine  32,  we  can  choose  any between  .007  and  .0035  (see  figure  9). 
Consequently,  choos ing  4^^  = .005,  we  have  that. 


— . ooSCrod) 

— .oo^h  = -.ooS'o  +-  <32 e 


(14) 


which  yields  an  82  = + .013.  Similarly,  for  a^  choosing  an  .0025, 

we  have  that. 


-.,0025(700)  -..OOiSC'S^s)  (15) 

ool  • — — • COS’D  4-  . 015  €(  + <35^ 


yielding  an  a^  = - .0093  °K/meter.  The  resultant  recovery  is  shown 
in  figure  10.  Except  for  the  location  of  the  inversion  points,  the 
recovery  is  nearly  exact.  In  this  method,  the  location  of  the  inversion 
points  are  critical  and  consequently,  enough  measurements  of  temperature 
with  must  be  made  to  locate  either  the  inflection  points  or  slope 
changes.  As  the  peak  of  the  oxygen  absorption  line  has  an  extinction 
coefficient  of  about  .003/meter,  and  = .0035ec  0,  then  a radiometer 
centered  at  the  peak  frequency  can  cover  a range  of  .021  to  .003  as 
it  moves  from  the  zenith  to  horizon  from  0°  to  82°.  For  values  less 
than  .003,  a change  of  frequency  would  be  necessary.  Assuming  a non- 
changing value  of  k^,  at©=  45°,  an  error  of  + 1°,  would  produce  an 
error  of  + .0001  in  4>e  = .0040. 


Consequently,  depending  on  the  value  of  CP^ , we  conclude  that  we  can 
differentiate  a /\(p^of  .0005  easily,  as  a result,  from  = .010  to 
.003  we  conclude  that  at  least  15  measurements  in  6 are  possible. 
However,  depending  on  the  profile,  15  measurements  are  not  necessary 
in  general. 

The  requirement  that  a radiometer  be  able  to  distinguish  temperature 
differences  can  be  determined  from  the  following.  In  determining  the 
apparent  a,  we  used, 

r To  + (16) 


If  we  assume  that  can  be  measured  accurately  and  that  <^>5 is  not 
limiting  then. 


Consequently,  the  tolerable  error  in  temperature  depends  on  the  4)o 
used  and  the  tolerable  error  on  a.  The  tolerable  error  on  a is  on  the 
order  of  5%.  Assuming  an  average  a of  .005  °K/meter  and  an  average <po 
of  .005/meter,  then  the  tolerable  error  on  the  temperature  measurement 
is  .05  °K. 

As  expressed  in  section  IV,  there  the  minimum  temperature  resolvable 
is  given  by. 


-r 


C CkJ-O 


(18) 


where,  the  symbols  have  the  previous  meanings.  To  achieve  a system 


60 


temperature  noise  of  .05°K  would  require  an  integration  time  of  10  sec. 
with  a bandwidth  of  60  Mflz  and  a system  noise  of  4 DB.  Although  this 
would  be  a sophisticated  system, it  is  within  the  state-of-the-art. 

In  order  to  determine  the  extent  of  the  degradation  of  the  profile 
recovery,  we  will  assume  an  error  on  the  slope  measurement  of  + 5%. 

In  addition,  we  will  assume  that  the  error  on  4^  is  not  limiting  and 
proceed  to  determine  the  recovery  possible.  It  must  be  noted  that  this 
is  a gradient  method  as  opposedto  a spatial  method  (e.g.  Lidar  back- 
scattering)  . If  there  are  no  inversions,  then  the  error  on  the  gradient 
is  just  f 5%  as  only  one  measurement  is  necessary.  For  example,  for  a 
typical  gradient  of  - .005  °K/meter,  the  error  at  2km  for  a base 
temperature  of  291.00  °K  is  just  + 0.5  °K.  If  there  is  one  inversion, 
then  the  location  of  the  inversion  position  is  critical  and  the 
error  is  greater.  For  the  profile  under  consideration,  we  impose  an 
error  of  + 5%  and  proceed  to  determine  the  extent  of  the  deterioration. 
At  an  a of  - .005  °K/meter,  a change  in  ^ of  + .0003  is  necessary  to 
distinguish  a real  difference  in  slope.  The  first  inversion  point 
(change  in  slope  of  figure  9)  could  occur  anywhere  between  a 4^^of 
.0075  to  .005  which  from  figure  8 corresponds  to  an  H of  650  and  950 
meters  respectively. 

As  a^  is  derived  from. 


Q ( oue  to  Qj)  r.  0 ; 


(191 


then , 


2:i  a 


- 4^0  H I 


H 


(20) 
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For  a - .005  and  an  = + .015  (actual  value)  the  error  in  ^ due 
to  an  error  in  the  placement  of  the  inversion  height  becomes. 


A a = 7. S'  X lo'^e 


where  we  have  used  an  error  of  + 150  meters  on  H.  Then  the  error  in 
the  apparent  gradient  is, 


- ± 2.0  K \ o 


This  error  is  consistent  with  our  previous  remarl<s  on  the  distin- 
guisbability  in  Similarly,  the  error  in  the  apparent  a due  to  a- 


Aa  -±  2.S  \ io’*<:5x)o-3)e 


(/oo) 


A a - i:  5.7  X I o' 


The  combined  error  for  a two  inversion  system  is  just 


ACL  - i 2.1  X I o' 


Consequently  , we  c-onclude  that  for  an  error  in  placement  of  the 
inversion  heights  of  + 150  meters  in  a two  inversion  system,  the  error 
in  ^ is  on  the  order  of  + .0002. 

From  the  analysis  above  we  can  now  construct  a degraded  profile. 

We  place  the  inversion  heights  at  the  extremes  (i.e.  first  inversion 
at  650  meters  and  the  second  inversion  at  1650) . For  the  first  gradient 
we  have  that, 


a 


Q|  -h  Q 2 ^ 


(27) 


then , 
or, 

^2  — . ooo  ^ ± . 0005^ 

(28) 


- 4-.  ot  O ii  . OC3  ^ 


(29) 


For  a^,  with  at  1650,  we  have  that, 
~.003(  ±,0002.  = -,ooS  ± .0002b-^j^OIO^ 


or. 


. OOl  q ±. 000'S  — ■i’.  0020  ±..  OOH>  -tOjC-Ol  (31) 


or. 
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(32) 


.ooo|  ±.  . oo  f (o)  Cl- 1 2 X to  0 - 


— — , oOt>2.  i'.OlO 


Using  these  values  we  constructed  the  degraded  profile  which  is  shown 
in  figure  10,  labeled  degraded  prof ile, using  the  average  values  derived 
above . 

Again  we  conclude  that  the  recovery  of  the  temperature  profile 
is  satisfactory.  Except  for  the  exact  placement  of  the  inversion  heights, 
the  recovered  profile  is  adequate  for  the  intended  use  in  radio-wave 
propagation . 

B.  Water-Vapor  Profile  /nalysis 

Although  originally  we  had  proposed  to  recover  the  water-vapor 
profile  spatially,  we  now  believe  that  the  application  of  the  gradient 
method  used  for  the  temperature  recovery  is  a more  powerful  technique 
requiring  less  measurements  and  leading  to  better  recovery. 

The  equivalent  water-vapor  density  integrated  along  the  line  of 
sight  is  given  by, 

VJ  - 

where  the  symbols  have  the  previous  meanings.  If  we  assume  that  P 

o 

and  Tq  are  the  conditions  at  h^  and  P and  T the  conditions  at  \\2  then. 


In  addition  if  we  linearize  the  density  profile  between  h^  and  h2  we 
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(3) 


yO  - +-  ^ Vi 


where  a is  the  gradient  between  h and  h . Substituting  we  have  that. 


(4) 


Integrating  we  have  that, 


/ P' 

1 ^ _ Pa  (■ 

( 

1 

(5) 


If  we  let  h 

1 


0 (i.e.  at  the  ground) then. 


(6) 


Consequently,  with  ground  conditions  known,  yO^  , P^,  Tq  and  the 

temperature  profile  known)  a measurement  of  V7  allows  the  determination 
of  the  gradient.  In  addition, as  the  expression  is  not  very  sensitive 
to  the  temperature  ratio  we  can  write, 


yg  Vo/p 


(7) 


The  ratio  P^/P  we  have  utilized  is  the  standard  pressure  gradient  as 
deviations  are  not  great. 

In  order  to  utilize  equation  (7) , we  have  choosen  increments  of 
200  meters  to  extract  10  measurements  to  probe  up  to  2 km.  Larger 
increments  could  be  used;  however,  due  to  the  fact  that  the  water- 
vapor  profile  tends  to  have  a greater  number  of  inversions  then  the 
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Assuming  a temperature  profile  is  known , (previous  example)  we  first 

calculate  the  saturated  water-vapor  profile  W,  consistent  with  the 

temperature  profile.  This  is  shown  in  figure  11. This  was  computed 

exactly  from  equation  (1).  Shown  also  in  figure  His  the  calculated  W 

for  the  water-vapor  distribution  given  in  figure  6. 

From  section  IV,  part  B,  figure  3,  we  can  calculate  the  ^ necessary 

2 

for  the  saturated  case  to  produce  a W of  .3  gms/cm  at  200  meters.  From 
figure  3,  ^ = 9.50.  For  the  actual  profile,  at  200  meters, 11.80 
for  W = .20  gm/cm^.  At  this  ^ , the  radiometer  would  produce  an  SNR 
consistent  with  a change  of  W=  .20  (figure  3,  section  IV,  part  B,  top 
curve)  . 

Proceeding  in  this  fashion  we  go  to  the  next  200  meters  and  so  on. 

Table  I shows  the  result  of  this  procedure. 

TABLE  I 


H 

(meters) 

W 

(pr . cms . ) 

4>. 

2000 

1.43 

3.95 

1800 

1.30 

4.15 

1600 

1.16 

4 . 40 

1400 

1.06 

4 .65 

1200 

.98 

4 .90 

1000 

.84 

5.30 

800 

.72 

5.75 

600 

.62 

6.20 

400 

.45 

7.50 

200 

.20 

11.80 

With  these  values  of  W,  we  can  reconstruct  the  profile.  For  the  first 
200  meters,  applying  equation  (7)  we  obtain  an  a = + .050  g/m^/  meter. 
Since^yOg  is  5 g/m^,  then  at  200  meters  it  is  15  g/m^.  Continuing  in 
this  fashion  we  obtain  Table  II.  The  resulant  profile  is  shown  in 
figure  12. As  can  be  seen,  the  recovery  is  remarkable.  The  only  major 

6G 


4 


Hei^t  H (meters) 


2000 


FIGURE  11.  CALCULATED  W FOR  SATURATED  AND  MODEL  PROFILES  (SEE  TEXT) 


6 


error  occurs  when  an  inversion  falls  between  a 200  meter  increment; 
however,  the  procedure  tends  to  average  and  recovers  the  profile  well 
above  the  inversion  point  (see  top  of  the  recovery  curve) . 

TABLE  II 


H 

(meters) 


0 to 
200 
400 
600 
800 
1000 
1200 
1400 
1600 
1800 


200 

to  400 
to  600 
to  800 
to  100 
to  120 
to  140 
to  160 
to  180 
to  200 


In  order  to  evaluate  the  expected  error, we  note  that  in  Table  I, 

we  have  tacitly  assumed  that  differences  on  the  order  of = 0.1 

can  be  measured.  However,  we  showed  in  section  IV,  Tablel,  that  the 

measurable  difference  in  ^ is  dependent  on  the  actual  profile.  For 

this  case  we  proceed  as  follows.  The  measured  W at  the  first  200  meters 
2 

is  0.20  g/m  . The  next  measurable  W at  an  SNR  = 1 can  be  determined 
by  consideration  of  the  following.  At  400  meters,  = .45  and  ^ = 
7.50.  Hence  0.23.  From  figure  3,  section  IV,  part  B , at  ^ = 

7.50  and  an— S=  0.23,  the  SNR  is  well  above  1 and  can  be  measured. 

At  600  meters,  ^2  ~ ^ ~ 6.20  with  — S = 0.12,  the  SNR  is  still 

above  1.  At  800  meters,  W2  = .72  at  5.75  with  — & = .07  which 

is  just  measurable.  Consequently,  above  800  meters,  the  recovered 
profile  will  begin  to  degrade.  The  next  measurable  W occurs  at  about 
1050  with  a W2  = .88  and  ^ = 5.20.  This  requires  a recalculation 
of  ^ . Continuing  in  this  manner  results  in  the  degraded  profile  shown 
in  figure  12  labeled  degraded.  As  can  be  seen,  except  for  the  locations 
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of  the  inversion  heights,  the  degraded  profile  follows  the  true 
profile  remarkably  well. 

It  is  to  be  noted  in  passing  that  the  gradient  technique  works 
well,  as  it  appears  that  it  can  tolerate  greater  errors  without 
degrading  the  true  profile  appreciably  more  so  than  a spatial 
technique. 
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VII  CONCLUSIONS  AND  RECOMMENDATIONS 

The  primary  motivation  for  this  study  was  to  determine  the 
feasibility  of  using  the  existing  technology  of  infrared  and 
microwave  radiometers  to  measure  the  verical  distributions  of 
atmospheric  temperature  and  water-vapor.  This  information  was 
necessary  to  correct  for  radio  wave  refractive  index  in  order 
to  adequately  predict  the  radio  wave  trajectory. 

This  study  has  shown  that: 

1)  high  spatial  resolutions  of  temperature  and  water-vapor 
distributions  are  not  necessary  for  predicting  the  radio 
wave  trajectory 

2)  for  angles  of  elevation  greater  than  30°  off  the  horizon 
with  a tolerable  error  on  angle  trajectory  of  .1°  or  greater, 
a standard  atmospheric  gradient  model  is  sufficient  for 
radio  wave  propagation 

3)  unless  intial  trajectories  of  the  radio  wave  are  close  to 

the  horizon  (less  than  1°) , the  vertical  distribution  of  tern-  j 

perature  and  water-vapor  obtainable  from  the  passive  probes 
are  more  than  adequate  for  radio  wave  propagation  corrections 

4)  temperature  profiles  are  obtainable  using  the  natural  micro- 
wave  radiation  from  atmospheric  oxygen  from  a gradient  tech- 
nique 

5)  water-vapor  profiles  are  obtainable  using  the  natural  rad- 
iation in  the  infrared  from  atmospheric  water-vapor  from  a 
gradient  technique 
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6)  spatial  resolutions  obtainable  by  both  methods  are  more 
than  adequate  to  describe  radio  wave  propagation  for  most 
applications 

7)  both  methods  are  capable  of  probing  for  temperature  and 
water-vapor  up  to  2 kilometers 

8)  low  visibility  does  not  interfer  with  either  probe  appre- 
ciably 

9)  suuficient  number  of  measurements  (i.e.  distinct)  can  be 
made  to  accurately  track  at  least  5 inversions  in  water-vapor 
and  at  least  2 inversions  in  temperature 

10)  no  new  technology  is  necessary  to  implement  the  fabri- 
cation of  a temperature/water-vapor  probe 

Although  the  mathematical  excursions  used  in  the  previous 
sections  to  determine  feasibility  have  been  rather  unwieldy,  the 
actual  implementation  of  the  technique  is  remarkably  straight 
toward.  The  necessary  computations  and  comparsions  that  must 
be  made  in  real  time  can  easily  be  accomplished  by  a dedicated 
microprocessor/controller.  This  system  must  be  interfaced  with 
the  radiometer  along  with  local  sensors  for  temperature,  pressure 
and  humidity.  If  properly  implemented,  the  gradient  technique 
appears  to  be  less  sensitive  to  experimental  errors  than  a spatial 
technique  either  active  or  passive  and  holds  the  potential  of 
becoming  a more  reponsive  probe  than  either  Lidar  or  passive  probes 
using  weigthing  functions. 

It  is  recommended  that  the  next  phase  of  study  should  be  a design 
phase  in  which  the  integration  of  a microprocessor/  controller  with 
local  sensors  and  existing  radiometer  be  examined.  As  a first  step, 
it  is  recommended  that  this  design  study  begin  with  an  existing 
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infrared  radiometer  to  be  modified  for  integration  of  sensors 
and  computer.  The  infrared  radiometer  should  be  examined  first 
as  it  is  believed  that  the  microwave  radiometer  system  may  re- 
quire more  of  an  effort.  Consequently,  we  recommend  the  follo- 
wing task  items  in  a future  study: 

1.  Begin  design  study  for  intecjration  of  pressure , temperature 
and  humidity  sensors  with  a microprocessor/contol ler . 

2.  Work  with  an  established  infrared  firm  (C‘.g.  Barnes  Engin- 
eering) to  custom  design  or  modify  and  existing  infrared 
radiometer . 

3.  Begin  design  of  microcomputer  (in  liaison  with  an  existing 
computer  firm  such  as  Mostek)  to  establish  dedicated  ROM's 
for  the  necessary  on-line  calulations  along  with  RAf-1 ' s for 
processinf  the  in-coming  data. 

4.  Integrate  efforts  with  infrared  and  computer  companies  to 
develop  a system  suitable  for  carrying  out  the  determination 
in  real  time  the  water-vapor  density  profile. 
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